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ABSTRACT 

Compton scattering of photons by nonrelativistic particles is thought to 
play an important role in forming the radiation spectrum of many astrophysical 
systems. Here we derive the time-dependent photon kinetic equation that 
describes spontaneous and induced Compton scattering as well as absorption 
and emission by static and moving media, the corresponding radiative transfer 
equation, and their zeroth and hrst moments, in both the system frame and 
in the frame comoving with the medium. We show that it is necessary to use 
the correct relativistic differential scattering cross section in order to obtain 
a photon kinetic equation that is correct to first order in e/rrie, Te/rrie, and 
V, where e is the photon energy, Te and rUe are the electron temperature and 
rest mass, and V is the electron bulk velocity in units of the speed of light. 

We also demonstrate that the terms in the radiative transfer equation that are 
second-order in V usually should be retained, because if the radiation energy 
density is sufficiently large compared to the radiation flux, the effects of bulk 
Comptonization described by the terms that are second-order in V are at least 
as important as the effects described by the terms that are hrst-order in U, 
even when V is small. The system- and fluid-frame equations that we derive are 
correct to hrst order in e/mg. Our system-frame equations, which are correct 
to second order in V, may be used when V is not too large. Our huid-frame 
equations, which are exact in V, may be used when U —> 1. Both sets of 
equations are valid for systems of arbitrary optical depth and can therefore be 
used in both the free-streaming and the diffusion regimes. We demonstrate that 
Comptonization by the electron bulk motion occurs whether or not the radiation 
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field is isotropic or the bulk flow converges and that it is more important than 
thermal Comptonization if > STe/rrie. 


Subject headings: plasmas - radiation mechanisms: thermal - radiation transfer 


1. INTRODUCTION 

Compton scattering of photons by nonrelativistic particles is thought to play an 
important role in many astrophysical settings, including the early universe (see Peebles 
1971; Sunyaev & Zeldovich 1980), clusters of galaxies (see Rephaeli 1995), active galactic 
nuclei (see Mushotzky, Done, & Pounds 1993), compact galactic X- and gamma-ray sources 
(see Pozdnyakov, Sobol, & Sunyaev 1983), and supernova remnants (see McCray 1993). In 
this process photons lose energy to the electrons or gain energy from their thermal and bulk 
motions, as a result of Compton recoil. 

Starting from the Boltzmann equation for photons, Kompaneets (1957) derived the 
so-called Kompaneets equation, which describes the time evolution of the photon energy 
distribution caused by scattering by thermal electrons when there is no bulk motion, the 
radiation held is perfectly isotropic, and the change in the energy of the photon in each 
scattering is small. The conditions assumed in deriving the Kompaneets equation are never 
strictly satished in astrophysical systems, since the radiation held is always anisotropic near 
their boundaries. In inhomogeneous systems, the radiation held may be anisotropic even in 
the interior. Also, in many astrophysical systems the scattering particles have substantial 
bulk motions. Thus, although the Kompaneets equation has been used extensively to 
treat astrophysical systems, in many cases it does not give accurate results. This has 
motivated several authors to derive photon kinetic equations that are valid under more 
general conditions. 

Comptonization in regions where the radiation held is anisotropic has been studied 
either by using Monte Carlo techniques or by solving partial diherential equations for the 
specihc intensity of the radiation held using hnite-diherence methods. Pozdnyakov et ah 
(1983) have reviewed Monte Carlo methods and results. Nagirner & Poutanen (1994) have 
reviewed work based on calculation of the complete Compton scattering kernel for polarized 
radiation. Babuel-Peyrissac & Rouvillois (1969), Pomraning (1973), Payne (1980), Madej 
(1989, 1991), and Titarchuk (1994) derived radiative transfer equations that describe 
Comptonization of an anisotropic radiation held when there is no bulk motion. Chan & 
Jones (1975), Blandford & Payne (1981a), and Fukue, Kato, & Matsumoto (1985) derived 
photon kinetic equations for thermal particles with non-zero bulk velocity, in the dihusion 
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approximation. Thorne (1981; see also Thorne, Flammang, & Zytkow 1981) derived 
fluid-frame moments of the radiative transfer equation in general relativity. The equations 
derived by these various sets of authors have been widely used to study Comptonization 
by strong shocks and accretion flows onto compact objects (see, e.g., Blandford & Payne 
1981b; Payne & Blandford 1981; Lyubarskij & Sunyaev 1982; Colpi 1988; Riffert 1988; 
Mastichiadis & Kylafls 1992; Miller & Lamb 1992; Titarchuk & Lyubarskij 1995; Turolla et 
ah 1997; Titarchuk, Mastichiadis, & Kylafls 1997). 

In the course of our investigation of the effects of Comptonization on the X-ray spectra 
of accreting neutron stars (see, e.g.. Lamb 1989; Miller & Lamb 1992; Psaltis, Lamb, & 
Miller 1995) we have rederived the radiative transfer equation and its moments for static 
and moving media and found important corrections to almost all the above derivations of 
the photon kinetic or radiative transfer equations, as we explain in § 2. There we show that 
it is necessary to use the correct relativistic differential scattering cross section in order 
to obtain a photon kinetic equation that is correct to first order in e/rrie, Te/rrie, and V, 
where e is the photon energy, and mg are the electron temperature and rest mass, and 
V is the electron bulk velocity in units of the speed of light (we use units in which the 
Boltzmann constant and the speed of light are equal to unity). In § 2 we also demonstrate 
that the terms in the radiative transfer equation that are second-order in V usually should 
be retained, because in many situations the second-order terms are at least as important as 
the first-order terms, even when V is small (see also Yin & Miller 1995). If the terms that 
are second-order in V are instead neglected, significant errors are introduced in the photon 
kinetic equation and its moments. 

In § 3 we state our assumptions and approximations and introduce our notation. 

In § 4 we derive the time-dependent photon kinetic equation that describes spontaneous 
and induced scattering by static and moving media, the corresponding radiative transfer 
equation, and their zeroth and first angular moments, in the system frame and in the 
frame comoving with the medium. We derive the moment equations as well as the kinetic 
and transfer equations because, although it is usually necessary to solve the full radiative 
transfer equation in order to determine accurately the angular distribution of the radiation 
held, the moment equations can be used to speed up the numerical computation by a large 
factor (Mihalas 1980; see also Mihalas 1978, pp. 157-158) and to provide additional physical 
insight. The system- and fluid-frame equations that we derive are correct to first order in 
e/TTig. Our system-frame equations, which are correct to second order in V, may be used 
when V is not too large. Our fluid-frame equations, which are exact in V, may be used 
when K —1. Both sets of equations are valid for systems of arbitrary optical depth and 
can therefore be used in both the free-streaming and the diffusion regimes. Our equations 
can easily be generalized to describe scattering by an arbitrary number of particle species. 
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In Appendix A we give the photon kinetic and radiative transfer equations that 
are obtained by averaging the equations for a single electron over a drifting, relativistic 
Maxwellian electron velocity distribution. In Appendix B we give the radiative transfer 
equation that describes absorption and emission in moving media, and its zeroth and 
hrst moments. There we point out that the addition of a photon source term in the 
transfer equation without any corresponding absorption term (see, e.g., Blandford Sz 
Payne 1981) is fundamentally inconsistent with thermodynamics and leads to a radiative 
transfer equation that has a different mathematical character from the thermodynamically 
consistent equation. 

Finally, in § 5 we summarize our results and their implications for Comptonization by 
static and moving media. 


2. MOTIVATION 

The radiative transfer equation that describes scattering of photons by particles is an 
integro-differential equation in which only derivatives with respect to the spatial coordinates 
appear (see, e.g., Nagirner & Poutanen 1994). The scattering kernel in this equation 
is non-local in photon energy and depends on the (possibly complicated) correlations 
between the angular dependence of the specihc intensity of the radiation held, the velocity 
distribution of the particles, and the differential scattering cross section. In order to 
accelerate numerical calculations, gain better physical insight, and facilitate comparison 
with previous studies that made similar approximations, we convert this integro-differential 
equation into a partial differential equation over the spatial coordinates and photon energy 
by expanding the scattering kernel in powers of the dimensionless quantities e/rue, Te/rrie, 
and V, which we assume are small compared to unity. In this way, the scattering kernel 
in the transfer equation becomes local in photon energy and the scattering integral over 
solid angle can be expressed in terms of the angular moments of the specihc intensity of the 
radiation held. 

We are primarily interested in deriving a transfer equation that can be used to calculate 
the spectra of X-ray and soft y-ray sources, so in expanding the scattering kernel we shall 
keep only terms of zeroth and hrst order in Aeq/eq ~ Ejvrif, <C 1, which is the average 
fractional decrease in the energy of a photon in a single scattering in the electron rest frame 
(see eq. below). Then, in order to obtain a radiative transfer equation of consistent 
accuracy, the terms in the expansion of the kernel in powers of Te/rue and V that are of 
the same size as the terms of hrst order in Aeq/eq, i.e., of the same size as E/rrie, must be 
retained. Depending on the situation, terms of diherent order in Tg/me and V may need to 
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be included. 

We now discuss the accuracy of the expansion of the differential scattering cross section 
needed to obtain a transfer equation that is accurate to hrst order in e/rrif.-, the orders of 
the terms in V that must be retained, and some subtle points that must be taken into 
account if the diffusion approximation is used. 


2.1. Approximate Scattering Cross Section 


In order to obtain a transfer equation that is consistently accurate to hrst order in 
e/me, it is necessary to use the correct relativistic expression for the differential scattering 
cross section in the frame in which one is working. We shall work in the electron rest frame 
and we therefore use the Klein-Nishina expression for the differential scattering cross section 
(see eq. |j^). To hrst order in e/me, the Klein-Nishina cross section is 


da 
dff Q 


Saj' 

Ibvr 


1+ ik ■ for 

1-2 —(l-/o-/') 


L 


( 1 ) 


where aT is the total cross section for Thomson scattering and Iq and Vq are the direction 
vectors of the photon in the electron rest frame, before and after the scattering. 

Use of the Thomson approximation for the diherential scattering cross section (Chan 
& Jones 1975; Payne 1980; Madej 1989, 1991) rather than expression (|^) introduces errors 
in the radiative transfer and moment equations of order e/mg, i.e., of the same size as 
the basic Compton ehect. Neglecting the angle dependence of the term of order e/me 
in expression (|I]) (Pomraning 1973; Blandford & Payne 1981a) introduces errors in the 
radiative transfer equation of this same order because, unlike the Thomson approximation, 
expression (|l|) is not forward-backward symmetric. Approximating expression (|I]) by the 
average of the cross section over the electron velocity distribution (Titarchuk 1994) neglects 
the effects of correlations between the angle dependence of the differential cross section, 
the specihc intensity, and the electron velocity, introducing errors in the radiative transfer 
equation of order e/mg, Tg/me, and V. Finally, use in the system frame of the Thomson 
or Klein-Nishina expressions for the differential scattering cross section (Pomraning 1973; 
Payne 1980; Madej 1989, 1991) introduces errors in the transfer equation of order V and 
Te/rrie, because these expressions are valid only in the electron rest frame. 
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2.2. Importance of the Terms of Order 

It has long been recognized that if the divergence of the electron bulk velocity is 
non-zero, the electron bulk motion causes a secular change in the energy of the photons 
(see, e.g., Chan & Jones 1975; Blandford & Payne 1981a; Fukue et ah 1985). The photon 
kinetic equations derived in these works do not include any terms of second order in V 
and predict that photons are systematically upscattered in energy by the electron bulk 
motion if the flow is converging. However, as we show with the examples that follow, 
photons are systematically upscattered by the electron bulk motion even if the terms that 
are first-order in V have, on average, no effect. Indeed, under some circumstances photons 
are systematically upscattered by the bulk motion even if the flow is diverging. 

To see this, consider a situation in which the bulk motion of the electrons can be 
described as isotropic turbulence in which the velocity correlation length is much smaller 
than the photon mean free path and the bulk velocity satisfies V ■ H = 0. If all the terms 
that are first-order in V are included correctly in the photon kinetic equation but terms 
that are second-order in V are not included (Fukue et al. 1985), the zeroth and first 
moments of the kinetic equation will indicate that the electron bulk motion has no effect 
on the photon energy distribution on time scales longer than the velocity correlation time. 
In reality, however, under the conditions stated the effect of the electron turbulent motions 
is completely analogous to the effect of electron thermal motions, that is, the turbulent 
motions cause the mean energy of the photons to increase steadily with time, at a rate 
proportional to the second moment of the turbulent velocity held. Hence, the electron bulk 
motion causes the mean photon energy to increase with time even though the flow is not 
converging; in fact, in this case no terms that are first-order in V have any effect on the 
radiation held, on average. If the turbulent velocity is high enough, the terms that are 
second-order in V will cause the mean photon energy to increase even in a diverging flow. 

Origin of terms second-order in V .—The origin of the terms in the photon kinetic 
equation that are second-order in V can be understood by considering a cold (Tg = 0) 
electron fluid with a uniform bulk velocity V in the system frame. For simplicity, let us 
analyze scattering of photons in the Thomson limit in the electron rest frame. In this 
limit, the angular distribution of the scattered photons is backward-forward symmetric in 
the electron rest frame, so the average photon energy in the system frame after scattering 
is 7 ^ {1 — {k ■ V)) times the energy before the scattering (see Rybicki & Lightman 1979, 
p. 198), where 7 = (1 — and the average is over the initial photon wave vectors k in 

the system frame. Consider first the case of a photon distribution that is isotropic in the 
system frame. In this case the average photon energy in the system frame after scattering is 
7 ^ times the energy before scattering, so in each scattering the average energy of a photon 
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is increased by an amount of order and hence there must be terms of order in the 
kinetic equation that describe this effect. In the more general case of a photon distribution 
that is anisotropic in the system frame these terms of order almost always produce a 
secular increase in the average energy of the photons, but there are also terms of order V 
in the kinetic equation that may cause the average energy of the photons to increase or 
decrease, depending on the velocity held and the angular distribution of the photons (there 
are generally terms of order V even if V ■ 17 = 0; see Fukue et al. 1985). 


Relative sizes of terms of different order in V .—The order in 17 of a given term in the 
various moment equations does not by itself determine whether the term should be retained 
in solving a given transport problem. This is because the terms of different order in V 
involve different angular moments of the radiation held, so the relative sizes of the moments 
must also be considered in determining which terms should be retained (see also Yin & 
Miller 1995). 


For example, in the zeroth moment of the radiative transfer equation, which is a scalar 
equation for the energy density of the radiation held, the terms that are hrst-order in V 
involve the scalar product of the vector quantity V with the vector quantity dehned by 
the radiation held, i.e., the radiation hux H (see eq. |^^). In contrast, the terms that are 
second-order in 17 do not involve H but instead involve the scalar quantity dehned by the 
radiation held, i.e., the radiation energy density J, and the radiation stress-energy tensor 
(again see eq. a: 


). In general, ViVjK'’^ is of the same size as 17^J. Hence, if J is 


sufficiently large compared to H, the terms in the zeroth moment of the radiative transfer 
equation that are second-order in 17 are generally at least as large as the terms that are 
hrst-order in 17. 


In the hrst moment of the radiative transfer equation, which is a vector equation for 
the radiation hux, the situation is reversed, in that the terms that are hrst-order in 17 
involve J and (see eq. [0) whereas the terms that are second-order in 17 involve H and 
Qpfc In general, VjVkQ^^^ is of the same size as 17(17 ■ H). Hence, one might be tempted 
to neglect the terms of order 17^ in the hrst moment of the radiative transfer equation if J 
is large compared to H, even if the terms of order 17^ are retained in the zeroth moment 
equation. However, this is in general unsafe, because it involves treating the same terms in 
the radiative transfer equation, from which the moment equations are derived, diherently, 
in taking the zeroth and hrst moment. In addressing a given transport problem, one can 
only determine which terms in the expansion in powers of 17 must be kept by considering 
the boundary conditions as well as the transfer equation. 

The radiation held quantities that appear in terms that involve only odd powers of 17 
are all of about the same size. Similarly, the radiation held quantities that appear in terms 



that involve only even powers of V are all of about the same size, although they may be 
much larger or smaller than the radiation held quantities that appear in the terms that 
involve odd powers of V. Hence, when H is -C 1, terms of order and higher may be 
safely neglected in the derivation of the transfer equation (see also Yin & Miller 1995). 


2.3. Use of the Diffusion Approximation 

In situations in which the change during a scattering of the energy of a photon as 
measured in the huid frame can be neglected and the longest photon mean free path Amax 
is much smaller than the smallest length scale on which physical variables change, the 
specihc intensity in the frame comoving with the electron fluid, //(s), is given to lowest 
order in Amax/hmin by (see Mihalas & Mihalas 1984, p. 457) 

I,{E)^Jf{E) + 3i-H,{E), (2) 

where Jf{s) and Hf{e) are the zeroth and hrst angular moments of !/{£)■ The diffusion 
approximation (sometimes called the zeroth-order diffusion approximation; see Mihalas 
& Mihalas 1984, § 97) consists in assuming that If{e) is given exactly by equation (H). 
Then the hrst huid-frame Eddington factor = Kj (e)/Jf{e) is exactly equal to T 

here Kj{e) is the second moment of As is evident from equation (^, if the source 

function involves only the zeroth, hrst, and second moments of //(£), then in the dihusion 
approximation it is only necessary to solve the three equations consisting of the zeroth 
and hrst moments of the huid-frame radiative transfer equation and the closure relation 
r/i^) = I for the three moments Jf{s), and Kj{e); the specihc intensity may then 

be calculated in this approximation using equation (^. 

The moments of the specihc intensity in the system frame do not satisfy equation i). 
even when Amax/Tmin -Cl. If expression (H) is mistakenly used to relate the specihc 
intensity to its zeroth and hrst moments in the system frame (Blandford & Payne 1981) 
rather than in the huid frame, then errors of order V will be introduced in the radiative 
transfer equation and its moments, even if Amax/Tmin -C 1 (see Fukue et al. 1985). 

If one wants to use the dihusion approximation (|^) to solve for the moments of the 
specihc intensity in the system frame, one must solve the four equations consisting of 
the zeroth and hrst moments of the system-frame radiative transfer equation and the 
two closure relations involving the hrst two Eddington factors for the hrst four moments 
J{s), H\e), and Q*-^^(e) of the specihc intensity I{e) in the system frame. The 

reason is that when equation (^ is boosted into the system frame and the terms that are 
second-order in V are retained, as in general they must be (see above), the third moment 
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of the specific intensity is introduced into the closure relation. (Fukue et ah [1985] 

were able to set up a closed system of equations consisting of the zeroth and hrst moments 
of the system-frame radiative transfer equation and a closure relation involving the hrst 
Eddington factor only because they neglected all the terms in the radiative transfer equation 
that are second-order and higher in V.) 

When, as in Compton scattering, the energy of a photon in the huid frame changes in 
a scattering, expression may not be accurate for all photon energies, even if the mean 
free path of a photon is independent of its energy and Amax/Emin -Cl. As an example, 
consider a slab that is inhnite along the x and y axes but hnite along the z axis, in which 
the electrons are cold and have no bulk motion. Suppose the slab is illuminated from one 
side with monochromatic photons of energy propagating in the z direction. Since the 
electrons are static, any photons that have scattered have lost energy to the electrons and 
hence no longer have energy Sin- Therefore all photons anywhere in the slab with energies 
equal to £in have never been scattered and are still propagating in the 2 ; direction. As a 
result, the angular distribution of these photons is not described accurately by equation (|^), 
even if Amax/Tmin is very small; indeed, the Eddington factor at energy Sin is equal to 
unity everywhere in the slab. In this example, expression (|]) is accurate only at energies 
sufficiently below the injection energy ein. 

The electron bulk velocity E is a vector quantity, and hence the relative sizes of the 
terms of hrst and higher order in V in the expansion of the transfer equation and its 
moments depend strongly on the relation between the angular dependence of the bulk 
velocity, the radiation held, and the diherential scattering cross section. Therefore, even 
when the dihusion approximation is only slightly inaccurate (as for example when the 
Eddington factor dihers only slightly from A), this inaccuracy produces additional terms 
in the moment equations that are of second order in V and that are therefore of the same 
magnitude as the second-order terms that would be present if the radiation held were given 
exactly by equation (^. 

For these reasons we derive the radiative transfer and moment equations without 
making use of the dihusion approximation. 


3. ASSUMPTIONS, DEFINITIONS, AND APPROXIMATIONS 

In the sections that follow we assume that the electron gas is nondegenerate (electron 
occupation number -C 1). For conciseness we shall refer to it as a huid, without 
implying anything about whether it is collisional or collisionless. We assume that photons 
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are scattered only by electrons, neglecting scattering by any other particles, and set 
h = c = ks = where h is Planck’s constant, c is the speed of light, and ks is Boltzmann’s 
constant. We indicate quantities evaluated in the rest frame of a particular electron by a 
subscript ‘0’ and quantities evaluated in the inertial frame momentarily comoving (locally) 
with the fluid, which we call the ‘fluid frame’, by a subscript ‘f’. The ‘system frame’ may 
be any global inertial frame (such as the frame at rest with respect to the center of mass of 
the system, if it is inertial). Quantities evaluated in the system frame have no subscript. 

Bulk velocity and temperature. —We dehne the fluid frame as the frame in which the 
energy and momentum flux density both vanish (Landau & Lifshitz 1987, p. 505). In this 
frame collisions between the electrons tend to establish an isotropic velocity distribution. 

The system-frame three-velocity -u of a given electron is related to its fluid-frame 
three-velocity -u/ by a Lorentz boost. The hrst and second moments of u in the system 
frame are 

(S) = V (3) 

and 

(u^) = (v^) + , (4) 

where V is the three-velocity of the fluid as measured in the system frame and 

V = u — V (5) 


is the peculiar velocity of an electron as measured in the system frame. 


If the electron momentum distribution in the fluid frame is a relativistic Maxwellian, 
then the second moment of the electron velocity distribution evaluated in the fluid frame is 


j J me 


mt 


( 6 ) 


and the second moment of the electron peculiar velocity distribution in the system frame is 

Te 


mp 


-W 


me 


( 7 ) 


Description of the radiation field. —In section 4.1 and Appendix A we derive the 
equation that describes the evolution of a radiation held interacting with a moving electron 
huid. We shall refer to this equation as the photon kinetic equation when it is written in 
terms of the photon mode occupation number and as the radiative transfer equation when 
it is written in terms of the specihc intensity of the radiation held; the two descriptions are 
equivalent (see, e.g., Mihalas 1978, p. 32), but the radiative transfer equation is more often 
used in astrophysical problems. 
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In deriving the photon kinetic eqnation we shall describe the radiation held by 
the nnmber n(l,e) of photons with energy e propagating in direction / with a given 
polarization state (we snppress the dependence on polarization state, because we consider 
only unpolarized radiation). The hrst few moments of n{l,e) are 


n = 

^ J n{i,e) dn , 

(8) 

n* = 

— / n{i,e)rdQ , 
dvr J 

(9) 

= 

^ J n{i, e) IT dn , 

(10) 

= 

— / nile)iHn^dn. 

dvr J 

(11) 


In dehnitions ®-( 0 ) the dependence of the moments on position and photon energy have 
been suppressed for brevity. Here and below we display the dependence of the photon 
occupation number on I and e in order to distinguish it clearly from its zeroth moment. 

In writing the radiative transfer equation we shall describe the radiation held by its 
specihc intensity = 2e^n{l,e)] here the factor of two accounts for the two photon 

polarization states. The hrst few moments of are 


J = 

An J 

f I (/, e) dn = 2e^n , 

(12) 

w = 

dvr J 

f I{i,e) rdn = 2£V , 

(13) 

= 

An J 

f I{i, e) IT dn = 2e^T^ , 

(14) 

= 

An J 

f I{i,e) iTl’^dn = 2£V^'^ , 

(15) 


where again we have suppressed the dependence of the moments on position and photon 
energy. 

Compton scattering .—In the rest frame of the electron the diherential cross section for 
scattering of unpolarized radiation is (see Berestetskh, Lifshitz, & Pitaevskh 1971, p. 297) 


da Sar f s'q ' 
dVtQ Ibvr \£o 


^0 ^0 


(16) 


where ax is the Thomson cross section, Sq and /q are the energy and direction of the incident 
photon, and 


= 


^0 


£0 

rUe 


(l-k- I'o) 


( 17 ) 
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and /q are the energy and direction of the scattered photon. 


The energy and direction of propagation of a photon in the electron rest frame are 
related to the same quantities in the system frame by the Lorentz transformations 


Eq = 7e(l — I ■ u) = 


7(1 + /o ■ u) 


and 


^"0 — — “i ^ + 

^0 


7 
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u 


(18) 


(19) 


The photon phase-space volume is Lorentz invariant (see Mihalas 1978, p. 495), i.e., 

ededVt = Eq dEo dflo , (20) 


so 


n(l,E) = no(lo,Eo) . 


( 21 ) 


Validity of the approximations .—In deriving the photon kinetic and radiative transfer 
equations and their moments in the system frame, we shall retain terms up to first order in 
e/me and second order in u (hrst order in Te/nie), neglecting terms of order {E/me)u or 
higher. These equations are therefore valid when 

— <1 and — (^2 + —(172 + —) <1. (22) 

me rOe \ nie J \ nie J 

These conditions are usually satished in accretion onto white dwarfs and neutron stars, 
but are not satished in accretion onto black holes, because the bulk velocity 17 —> 1 at the 
horizon. 


Our expressions for the photon kinetic and radiative transfer equations in the huid 
frame are correct to all orders in 17 but only to hrst order in Te/mg and Ef/rUe. These 
equations are therefore valid when 


— 1 and 

rrie 


and can be used where 17 is ~ 1. 


rUe 




< 1 


( 23 ) 


4. PHOTON KINETIC AND RADIATIVE TRANSFER EQUATIONS 
AND THEIR ZEROTH AND FIRST MOMENTS 
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4.1. Photon Kinetic Equation 


The photon kinetic equation in the system frame is 
'• d^k' 


k^d^n{k) = 


(2i)s 


d^p {k'p^ kp'^ n(k') 1 + n(k) f{p) 1 — f{p') 

—W {kp, k'p'^ n(k) 1 + n(k') f{p) 1 — f{p') | , (24) 


where k^ = {e,el), n{k) is the photon occupation number, p is the electron momentum, f{p) 
is the electron momentum distribution, W (^kp, fc'p') is the transition rate for the scattering 
k + p ^ k' + p', and in writing the total derivative on the left side we have used the Einstein 


summation convention. The left side of equation (|2^) is manifestly covariant. The ride side 
is also covariant, if the appropriate transition rates are used. 

It is convenient to integrate hrst over the photon states and then over the electron 
states, because the angular integrals are then much simpler. We will therefore assume for 
the moment that all the electrons have the same momentum p. Because each side of the 
photon kinetic equation (1^) is covariant, the left side can be evaluated in the system frame 
and the right side in the electron rest frame, at the wave vector fco that corresponds to 
the wave vector k of the photon in the system frame. The resulting equation for electrons 
moving with velocity u = p/{'ymf.) is (compare with Peebles 1971, p. 204) 


[dt + rd^ en{i, s) = SoUeo 


’^o(^l)Wo) 1 +^o) 


- Mk, £o) + no(/o, 4) 



(25) 


where dt = d/dt, di = d/dx\ the x* are the spatial coordinates, Ueo is the electron density 

and 


in the electron rest frame, Iq and Eq are related to I and £ by equations (|T^) and 
£q is related to Eq by equation (pT[) . The factor preceding the cross section in the hrst term 
of the collision integral is the Jacobian that corrects for the different phase spaces of the 
incident and scattered photons. 


We evaluate the collision integral in equation (^^ by hrst using the Lorentz invariance 
of the photon occupation number to relate ?7,o(/o,£o) to n{l\£') and relating n{i',£') to 
n(l',£) by expanding n(l',£') to second order in u and to hrst order in Elme-, which gives 
(see Peebles 1971, p. 204) 


^o(C4) - n{i\£) + 


{J'o ~ Iq] ■ u + [Iq ■ if] (Iq — I'q] ■ u H- (l — Iq ■ Iq 

Tflg 

\ [{^0 - ^ o ) ■ «] ^ > 


Ede n(l' ,e) 


(26) 
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where = d/de and d/ = d’^ jde^. We then nse this result, the Lorentz invariance of the 
photon distribution, and {da/dfl')o and Ueo^o expanded to second order in u to obtain the 
approximate photon kinetic equation for electrons moving with velocity u 


(dt + n(/, e) 


Sn^aT 

Ibvr 


I dfl' [l + (fo ■ I'of] [l - 2;^ (l - ^0 ■ /o) 
1 - (Jo ■ uj + (Jo ■ m) -U^ 

(Jq -io) ■u + 2 (Jo ■ uj (Jo - ^o) ■ U 


ed. 


2 L 


(Jo ~ io) ' u e^d/ H-fl — /q ■ I'o) (4 + sd^) 

Tfie 

+2— (l - io ■ I'Jj n{i, e) (2 + ede)\n{l', e) 

TTlfi ) 

\ - 2—) (l - / ■ n) n(l, e) , (27) 


-UeO-T 


which describes the effects of scattering by electrons with velocity u. 


In Appendix A we give the photon kinetic and radiative transfer equations that are 
obtained by averaging equation (|^) over a drifting, relativistic Maxwellian electron velocity 
distribution (see eqs. [@]-|0). The moment equations derived in the next two subsections 
can be obtained by computing the zeroth and hrst moments of the equations given in 
Appendix A. Here we follow the simpler approach of hrst computing the moments of 
equation (|27|) and then averaging them over the electron velocity distribution. 


4.2. Zeroth Moment and Radiation Energy Density 


We compute the zeroth moment of the photon kinetic equation by hrst integrating 
both sides of equation (p7|) over all directions. Making use of the Jacobian 


we hnd 



/ So deo\ 

/s' de'\ 





d{no,n') 

\e de ) 

\e'ode'o) 






1 - 2 (Jk 

■ v}j -\-2 

J' ■ u 

) -2u^ + 3(1 






+3 ^ 

i'-uf-A (Jo 

■ u'j (J' ■ uj , 

(28) 

3(Tt" f 




r 2e / 




1 + (^0 ■ 

1 - (1 - 

L rUe ^ 

k ■ k) ■ 



dtu + djU* 
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0 / ■^ 

\ 2 /- 

. \ 



/ ^ 

. \ 1 

[{ 

1 - 3m2 + 3 (/ 

’■u) +2{l 

'■“)] 

+ (/ 

o-«) 

1 + 2 

'■“)] 


edp 


— • m) I 3 + 6 (J' ■ uj + 1 + 2 (J' ■ ede + {i'q ■ uj e^dl | 

+ (^"o ■ "w) ^6 + Aed^ + ~ 


+ — (l-l ■ /') n{i, e) (2 + 

TTT/g \ / 


-rieCTr 


1 — 2 - ) 77, — WUi 

rrif 


n{i', e) 

(29) 


Next we integrate over dflo and then over dQ', using definitions (^-(pd]). The result is 


Hear 


n 


[ptn + dird^ = (3 + ede) n^Ui + (4 + + y {Aede + £^5^) 

/36 34 ^ 11 2 ^ 2 \ fa 1 2 

+ ned^n — n’‘ed^n^ + rd^ed^ind^ — , (30) 

where we have again used the Einstein summation convention. Finally, after averaging 
equation (^0|) over the electron velocity distribution we obtain the zeroth moment of the 
kinetic equation that describes the effects of scattering by a fluid of electrons, namely, 

{dtn + cljn*) = (3 + ede) n'Vi + {A + ede) n 


UedT 


rrie 


+- ((v^) + (Aede + e^dl') n 

+ (— + —+ —e'^df] {rd'^{viVj) — 

VIO 10 20 V V ^ 3 ' V 

+ ("^ + ^ede + — (n^WiVj - -nvA 
VIO 10 20 V V ^3 J 

+ - f—) ( 2^2 - 2n*n* + 

2 Vm-e/ ^ 

+ nedeU — n^edeU^ + ■nd^eden''^ — rd^^ederd^^^ . (31) 


This equation is valid in both the diffusion and free-streaming regimes, for any arbitrary 
(possibly anisotropic) distribution of electron velocities. 
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If the electron velocity distribution in the fluid frame is a relativistic Maxwellian with 
temperature T^, then (see eqs. [|| and [0) 


(viVj) ~ ~ —6ij , 

O TTl/p 


(32) 


and the zeroth moment of the photon kinetic equation can be written, to the same accuracy 
as equation (^), as 


(dtn + diTi"^ = (3 + ede) n^Vi + ^ (4 + ede) + (^ede + 


UeCTT 


/36 34 


£ 
me 
11 


n 


1 — + ) ( ri^^ViVj — -nV 


VlO 10 


20 


1 


4{^){2n^-2n‘n‘ + 2n‘V^-2n‘^^n‘^^ 

+ nedeU — n^edeTi^ + n^^edeTi^^ — , 


(33) 


where we have used the relation = n. The corresponding zeroth moment of the 

radiative transfer equation with emission and absorption included (see Appendix B) is 


dJ d 
dt dx^ 


{ede 

r/e-4Te\ 1 

- 

1 

L\ rUe / J 




H- Ede{6j) + sdeH^Vi 

rUe 


-AsdeJ + edl{eJ) 


+ ( —ede + —e^dt ) ( K^^V.Vj - -JV 


10 


20 


4 Vm. 


{edeJ - J) 4 - ^ 






+ [edeK^^ - — - {edeQ^^^ - 


UeUT 


+ (i - + irvs;) 


- ■! IJ - H‘V, I + 


K^^ViV. + -JV^- H% 


1 


ede 


+-K^^V,V^e^di\xe. 


(34) 


where rje and Xe are the emission and absorption coefficients, which are defined in the ffuid 
frame but evaluated at the energy £ of the photons in the system frame. 
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4.3. First Moment and Radiation Flux 


We compute the first moment of the photon kinetic equation by multiplying both sides 
of equation (^) by I and then integrating over all directions. In performing the integration 
we use the transformation (W) and the Jacobian d^), integrating hrst over cIQq and then 
over dQ'. The result is 


1 


rieCXT 


2 £• 1 

(dt'n’'+ = —n* — -—{1 + ed^) n’'+ - (^—nui + 3n’‘^Uj^ 


— - (Sn^^’^UjUk — 2n^UjUjj — — 3ede (nui) + ed^ (n^^Uj^ 


1 

To L 
1 

To L' 

3 


lOede (n^Uj^ Ui + Oedg — eds'ndu^ 

3e^dl (n^Uj^ Ui + e^dl 


+ ndedeU - . (35) 


Averaging equation (|55[ ) over the electron velocity distribution, we finally obtain 

—^— (dtn^ + djU^A = —n* —(1 + edg) n* 

rieCXT ^ ^ orrie 

+1 (n'’Vj - InrJ - (n‘W, + 3nVi') 

(4n«Hn4 - n'V,\4) 

-TeS. [9n‘^'‘VjVt + 10iiW,l4 - n‘VjVj) 

(n«‘i4Vi + srfV'v;) 

(An^^’^ivjVk) -n^{vjVj)^ 

(n'^^ivjVk) + lOn^VjVi) - n^VjVj)'^ 

(n^^’^ivjVk) + 3n^{vjVi)^ 

+ ^ (^) (2W - 

+ n^edsU — rC^eden^ + . (36) 


If the electron velocity distribution in the fluid frame is a relativistic Maxwellian with 
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temperature Te, then 


and 


nHvjVi) = —n* 
me 

(37) 

"^'"{vjVk) = —n* , 
me 

(38) 


and the first moment of the photon kinetic equation reduces to 
1 


UeO-T 


(dfU^ + = — 


1 + (1 + ^e) + (i + 

1 


5 rup 


n 


+- [n^^Vj - - —ede (n^^Vj + 3nld) 

-H (4n«'=l/,V4 - n‘VjV,) 

-—eds (9n^^%Vk + lOn^VjVi - nWjVj) 

+ rCeden — n^^ede'n? + . (39) 


The corresponding first moment of the radiative transfer equation with absorption and 
emission included is (again see Appendix B) 


dtH^ + d.K^^ = 


]ji _ 2 {Tp 3e 


rrip 




5 L 


edp 


e-m 


nip 


-' + I^ed^^e 

me 




1 


10 L 
1 


(9 - ede) K^^Vj + (7 - 3ede) JVi 


+17; (-6 + Sed, - E^df) WV,V. + - (1 + ec),) H-V,V, 


10 

1 


^\{^) bv-+)f-04//-i/0 




+ 

1 , 


(edeK^^ - K^’^) ^ - {edeQ^’^^ - 


Rijkl 


>ne(JT 




(l + W - KW, + ede 


e I Xe + (2 - ede) Vs , 


(40) 


where is the fourth moment of the specihc intensity. 
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4.4. Equations in the Fluid Frame 


The photon kinetic and radiative transfer equations as well as their moments take their 
simplest form in the fluid frame, since 1/ = 0 in this frame, by dehnition. These equations 
can be written choosing as independent variables either Eulerian coordinates hxed in space 
and time or Lagrangian coordinates comoving with the fluid. 


In terms of the fluid-frame coordinates, the transfer equation 
Erst moment (|34f) and (HO) in the fluid frame become 




and its zeroth and 


{dtf + l}dif)lf{if,ef) = Uejar 


-j (C2Jf + C^H) + CiKj + Ci'^Qf + 


{J, - i)H) + fflKf - 


(41) 


dtfJf + di^Hj — nej(JT< 


and 


Q I ~ 4Te\ Tg 2 


me. 


me 


J 


f 


iff/ 

4 \me 


{e,amf - i) 4 - - flf) ^ 

ijk' 

f 


+ {e,d,,K‘i - Kf) - {e,d„Qf - Qf) ^ 


-f 


f 


d,^m + d,^Kj = -nejaT\H} + - 


5 L 


^ /Ef 4Te\ Tg 2 ^ 


>(42) 


HI 


4 \me 


{e,a„H} - Hj) H 

C’ f 


f 


£ 


f 




/^ijk -nijkl 

Jk _J^ 


f I e) 


-f 


(43) 


where the coefficients Ci-Cq, which are given in Appendix A, are to be evaluated at F = 0. 

The left sides of the above equations can also be written in terms of the system-frame 
Eulerian coordinates as in equations (95.9), (95.11), and (95.12) of Mihalas & Mihalas 
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(1984). The resulting equations, which are often called mixed-frame equations, are correct to 
all orders in V and therefore can be used in situations where the bulk velocity is relativistic. 
The velocity-dependent terms that appear in the mixed-frame equations arise from Lorentz 
transformation of the left sides of equations (E3)-(E3), whereas the velocity-dependent 
terms that appear in the system-frame equations arise from Lorentz transformation of the 
scattering integral on the right sides of these equations. 

The right side of the radiative transfer equation can also be used in the formalism 
developed by Thorne (1981), for solving Comptonization problems in general relativity 
(compare eqs. [6.10] and [6.13] of Thorne 1981) 


5. DISCUSSION 

In the previous section we derived the radiative transfer equation in both the system 
and fluid frames, taking into account absorption and emission as well as spontaneous and 
induced Compton scattering. In this section, we first show that our equation reduces to 
the Kompaneets equation in the appropriate limits and call attention to several errors 
and misunderstandings in the literature. Next, we discuss the moment equations for an 
anisotropic radiation held in a static medium and show that the radiation force, and hence 
the critical hux and luminosity, generally depend on both the photon energy spectrum and 
the electron temperature. We then consider the equation for the zeroth moment of the 
specihc intensity for moving media and show that if the radiation held is isotropic, the 
terms in the transfer equation that are second-order in the electron bulk velocity produce 
a systematic increase in the energy of the photons that is completely analogous to the 
systematic increase in the energy of the photons produced by the electron thermal motions. 
We also show that Comptonization by electron bulk motion occurs, whether or not the 
radiation held is isotropic or the bulk how converges, and give estimates for the time scales 
on which the photon energy distribution changes because of systematic downscattering and 
upscattering caused by the electron thermal and bulk motions. We derive a new, more 
general condition for determining when Comptonization by the electron bulk motion is 
more important than Comptonization by the electron thermal motions. We conclude by 
indicating how the transfer equations we have derived can be solved using the method of 
variable Eddington factors. 
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5.1. The Kompaneets Equation 


When the radiation field is isotropic and there is no bulk velocity and no absorption or 
emission, equation (^) reduces to the Kompaneets (1957) equation^ 


dtj = 


F AT T 

— ds{eJ) - -ed^J H— -edl{eJ) + 

me me me 


— {edeJ - J) 

me 



(44) 


where on the left side we have introduced the differential Compton time dtc = ne<JTdt. 
The first two terms on the right side of equation (|^) describe the effects of systematic 
downscattering and upscattering of the photons by electrons. The third term describes 
the diffusion in energy produced by the thermal motion of the electrons. The last term 
describes the effect of induced Compton scattering. 


As noted by Kompaneets (1957), any Bose-Einstein distribution is a stationary solution 
of equation (|^). Induced Compton scattering cannot change the chemical potential, 
because it does not change the number of photons. Hence, the statement by Pomraning 
(1973, p. 193; see also Rybicki & Lightman 1979, p. 209) that the Planck spectrum (the 
particular Bose-Einstein distribution with zero chemical potential) is the only stationary 
solution of the Kompaneets equation if induced scattering is included is not correct. 
Integrating equation (^) over energy shows that if the electron temperature is equal to the 
Compton temperature 


Tc = 


- jne^) 

4{e) 


(45) 


the photon energy density remains constant although the photon spectrum may evolve with 
time; here the average is over photon energy, using the photon-number density N{e) = {J/e) 
as the weighting function. Note that, for a given photon spectrum, induced Compton 
scattering always decreases the Compton temperature. 


^ In their derivations of the Kompaneets equation, Rybicki & Lightman (1979, p. 213) and Katz 
(1987, pp. 100-114) did not take into account in the collision integral the different phase spaces of the 
incident and scattered photons. However, following Kompaneets they evaluated the integral by using 
photon conservation and the thermodynamic equilibrium photon distribution rather than by performing 
the integration directly and thereby obtained the correct result despite this error. These authors also used 
the Thomson approximation to the Klein-Nishina cross section. In general this introduces an error of the 
same size as the systematic downscattering term (see below), but this error vanishes if the photon distribution 
is isotropic. 
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5.2. Implications of the Moment Equations for Static Media 


As we mentioned in the Introduction, the Kompaneets equation (Q) is not strictly 
valid for astrophysical systems, since it requires the radiation held to be isotropic and hence 
that no radiation leaves the system. If the radiation held is anisotropic but there is no bulk 
motion, the system of equations 


and 


dtJ + d,W = v-xJ 

+neaT 


ed. 


e-m 


me. 


+ —edle 

me 


J 


(46) 


dtW + d.K^^ = -nearH^-xH^ 


--rieO-T 


Te - 3£ 


me 


+ ede 


ATe 


me 


+ ^edle 

me 




(47) 


must be solved simultaneously together with closure relations for and J, which are 
usually introduced as variable Eddington factors (see Mihalas 1978, p. 157). For simplicity, 
we have neglected the ehects of induced Compton scattering in equations (^61) and ( 0 )- 
The term {2ne<JT/^)[{Te — 3e)/me]i7* does not appear in the hrst moment of the transfer 
equation derived by Chan & Jones (1975), Payne (1980), and Madej (1989, 1991), because 
these authors did not use the appropriate relativistic scattering cross section. 


In a static medium and in the absence of absorption, emission, and induced scattering, 
F*, the radiation force per unit volume on the electrons (see also Miller & Lamb 1995) can 
be obtained by integrating equation (|47|) over photon energy, which gives 




dvr / 8 




me 




(48) 


where {s)^ is the average photon energy, using 77® as the weighting function. Blandford & 
Payne (1981a) neglected the terms in the square brackets on the right side of equation (^7]) 
and therefore the radiation force given by their equation is incorrect to hrst order in ejme 
and Te/me- Equation ( ^8]) shows that the volume radiation force on the electrons produced 
by Compton scattering of the photons is different from that obtained in the Thomson limit, 
and depends both on the photon spectrum and on the electron temperature (see also Fukue 
et ah 1985). Hence, to hrst order in s/me and Te/rUe, the energy-integrated critical radial 
radiation hux that produces a radially outward radiation force which exactly balances 
the inward gravitational force of a massive object depends on the photon spectrum and the 
electron temperature. For a completely ionized hydrogen gas, the critical radiation hux at 
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radius r is given implicitly by the equation 




crit 


crUpCM A _ 8 2 '^^ 

r'^ar I 5 me rrie 


(49) 


where rup is the proton mass, M is the mass of the object, and {e)jrcTit is the average photon 
energy, using as the weighting function. For example, in the spectral formation region 
of many neutron-star low-mass X-ray binaries, {e) is ~ 1 keV and Te can be ~ 25 keV 
(Miller & Lamb 1992; Psaltis et ah 1995), in which case the critical flux is ~ 10% smaller 
than the usual Eddington critical flux computed assuming Thomson scattering. Note also 
that to this order the critical luminosity Lent = generally depends on radius, 

because both (e) jrcrit and Tg generally depend on radius. 


5.3. Implications of the Moment Equations for Moving Media 


Consider now a medium in which the electron bulk velocity is not zero. Suppose 
hrst, for simplicity, that the radiation held is isotropic in the system frame. To make 
this situation concrete, consider a thought experiment in which electrons are moving with 
uniform and constant bulk velocity V through a box with sides of length L. The electrons 
are assumed to be able to pass through the walls of the box whereas photons are conhned 
inside the box and have a mean free path much larger than L. Under these conditions, 
scattering of photons by the walls of the box, which is much more frequent than scattering 
of photons by the electrons, keeps the photon distribution nearly isotropic. 

In this situation the equation for the zeroth moment of the specihe intensity reduces to 

+ edle^ J , (50) 


where for simplicity we have neglected absorption, emission, and induced scattering. 
Equation (|5(]| ) shows that when the radiation held is isotropic in the system frame, 
Comptonization by the bulk motion of the electrons is described entirely by terms that are 
second-order in V; all terms that are hrst-order in V vanish identically. 


Suppose now that (i) photons with energies e Tg + are injected into the box 

with an isotropic momentum distribution, (ii) the photons are allowed to remain in the box 
for a distribution of residence times that decreases exponentially with the residence time, 
and (hi) 


T, 


rrig 




yb = 4 


Les ^ 1 5 


(51) 
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where tres is the mean residence time measured in units of the Compton time (l/rigaT)- 
The solution of equation (pOD is then a power-law spectrum with a high-energy cutoff. At 
high and low energies the spectrum is 

t , £ -C Te + \meV‘^ , 


J = 


exp[—e/(Te + , e 3> Te + ImeV^ 


(52) 


where 


a = — 


(53) 


5-A + l. 

and the factor arises because we are considering the energy density rather than the 
photon occupation number. This solution is a simple generalization of the solution to the 
Kompaneets equation obtained by Shapiro, Lightman & Eardly (1976; see also Rybicki 
& Lightman 1979, Ch. 7) and shows that the terms in the transfer equation that are 
second-order in the electron bulk velocity produce a systematic increase in the energy of the 
photons that is completely analogous to the systematic increase in the energy of the photons 
produced by the electron thermal motions. In treating astrophysical systems, the mean 
residence time t^es in equation (|5T|) is often expressed in terms of r^, where r is the electron 
scattering optical depth of the system. For example, if photons are injected at the center 
of a uniform-density spherical cloud of optical depth r, then the mean residence time is 
(3/7r^)(r -|- 2/3)^ (Sunyaev and Titarchuk 1980). 

When the photon momentum distribution is not perfectly isotropic, the characteristic 
time scales on which the photon distribution changes because of systematic downscattering 
and Comptonization by the electron thermal and bulks motions are (see eq. [Q) 


t. 


-1 

down 

-1 
th 
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tv 
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-1 

V2 


UeCTT 


ne<JT 


HearV 


UeUT 



(54) 

(55) 

(56) 

(57) 


where the last two time scales have been estimated from the terms that are hrst- and 
second-order in V. 

Comparison of rates (^) and (^71) shows that Comptonization by the electron bulk 
motion is more important than Comptonization by the electron thermal motion if 


( 58 ) 
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In fact, Comptonization by the bulk motion occurs whether or not the radiation held is 
isotropic or the bulk motion converges (see also eq. [Q), contrary to the impression given 
by Blandford & Payne (1981a). Comparison of rates ( |5^) and (|F7|) shows that if J is 
sufficiently large compared to H^, the effects of bulk Comptonization described by the terms 
that are second-order in V are dominant compared to the effects described by the terms that 
are first-order in V. 


In estimating the characteristic time scale for bulk Comptonization, Blandford & Payne 
(1981a; see also Blandford & Payne 1981b and Payne & Blandford 1981) used only one of 
the several terms in their equation that are hrst-order in V, neglecting other terms that 
generally also produce systematic upscattering or downscattering of photons. As a result, 
the characteristic time scale that they obtained for bulk Comptonization is proportional 
to (V ■ C)/3, which they assumed to be proportional to n^aTV"^/?>. This assumption is 
not generally valid. When it is, their expression suggests that bulk Comptonization is 
less important than thermal Comptonization if < 12Te/me. However, comparison of 
rates (l5^) -(|57| ) shows that the bulk Comptonization terms that are second-order in V are 
already as important as the thermal Comptonization terms when ~ ST^/me and hence 
that these terms can be more important than the thermal terms even if < 12Te/me. 


5.4. A Method of Solving the Transfer Eqnation for Moving Media 

In more realistic models, the radiation held is not isotropic in the system frame in the 
presence of bulk motion (see, e.g.. Miller Sz Lamb 1993, 1996) because electron scattering 
tends to isotropize the photon distribution in the fluid frame. In this case, the full radiative 
transfer equation must be solved. This can be done using the method of variable Eddington 
factors (Mihalas 1980; see also Mihalas 1978, pp. 201-203), in which the radiative transfer 
equation and its zeroth and hrst moments are solved iteratively. In this approach, the 
second and higher moments of the specihc intensity are related to the zeroth and hrst 
moments via variable Eddington factors. The zeroth and hrst moments of the radiative 
transfer equation are then solved using initial guesses for the variable Eddington factors 
and the source function is computed from the calculated moments of the specihc intensity. 
The radiative transfer equation is then solved, the Eddington factors are updated, and the 
whole procedure is repeated until convergence is achieved. 

A detailed study of the solutions of the equations derived here for realistic models of 
astrophysical systems will be reported elsewhere. 
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APPENDICES 


A. PHOTON KINETIC AND RADIATIVE TRANSFER EQUATIONS 
FOR SCATTERING BY AN ELECTRON FLUID 


In this appendix we give the photon kinetic equation and the corresponding radiative 
transfer equation for scattering by an electron gas with temperature and bulk velocity V. 


We start from the photon kinetic equation (|27|) , transform all the quantities to the 
system frame using equations (p!^)- ([I^) , and then average over the electron velocity 
distribution. The resulting photon kinetic equation is 

^ (dt + Pdi) n(l, s) = e) + j (7^2r^ + 
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+]-edl\2^lH^+ v,v, + (i-vf rA , (A5) 
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2 L me ^ ^ J 

= lOrPVkVi + ^ed.rPVkVi + ^e^d^^JH^VkVi. (A7) 

The corresponding radiative transfer equation is 
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4tc \ TTZg / 

where 

= —1 + 2 -h f/ ■ , (-^9) 

TTlg 

C 2 = l- — + —ede + 2— + 2(i-V)+4(i-vY-V^ 

rrie rrie me ^ ^ ^ ^ 

■ (' ■ - 2 (' ■ ^d1+r'®' [ 2 ^+('' ^)1 • 

. / / tg J ^ L '' 

A = —r - —ed,p - 4—r -17 - 2 f/ ■ f) r - 8 (i ■ fr + 21 ^^/* 

me me me ^ ^ ^ ^ 

+ede 2—P + Vi + 2(i ■V)Vi + 2 (l -vfp 

- me 

+ l-£d^\-2^P -2(l-v)v] , (All) 

2 L me ^ 2 ' J 

Cl = PP - —PP + —edePP - Q—PP + 4 f/ ■ V) PP - 2PV. 

me me me ^ 2 

+10 (/ ■ vfpp - 8 (/ ■ f) rv^- - 3V‘^PP + ViVj 
+ede \-2—PP - (l-V) PP - 4 (l-vYpP 

l me ^2 V / 

+-edf h^PP + ViV. + (I-vYPP] , 

2 L me ^ V y J 

Cf = ^PPlk _ ^^Qjipik 

me me me 


(A12) 



pijkl 

^6 


+rF14 - ‘2l%Vk + 4 (f ■ v) iH^Vk 


+ede 

+-e^dl 
2 ^ 




rPi^ + rPVk + 2 (/ ■ - 2rVjVk 

-2—PPP-2(i-V) PPVk 

m.. \ J ^ 


PPVkVi + 2edePPVkVi + -e^diPVVkVi . 

2 


(A13) 

(AM) 


B. EMISSION AND ABSORPTION PROCESSES 
IN THE SYSTEM FRAME 

In the absence of scattering, the transfer equation in the system frame becomes simply 
(see, e.g., Mihalas & Mihalas 1984, p. 422) 


(dt + Pd,)l{i,e) = h(£/)- (j) 


(Bl) 


where again the subscript ‘f’ refers to quantities in the fluid frame and we have for simplicity 
assumed that the emission coefficient r] and the absorption coefficient x are isotropic in this 
frame. In writing equation (|Bl| ) we have also used the fact that the absorption and emission 
coefficients are dehned in the fluid frame and that {Tj/e^) and {ex) are Lorentz invariants. 

We can now expand ri{ef) to second order in V as 
T]{ef) = r][e+{ef-e)] 


- he + (£/ - £)deVe + 2^^ f - ^ehe 

^ he + -l-V^ edeTje + (/' v'f e^dlrje , 


(B2) 


where = ri{e). We also expand x{^f) ia a similar way. The transfer equation to second 
order in V then becomes 


{dt + l%)l{l,e) = {[l + 2{i-V) + 3{l-V)^ -V^ 


+ 


■{l-V)-2{l-V)^+ -V^ 


ede + -{i-Vfe'^dl \ r]^ 


l-{l.V) + -V^ 


+ 


-(/■!/) + (/-17)2 +-1/2 


edr 


+-{i-V)^eds}xeI{l,e). 


(B3) 
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We can then obtain the zeroth and hrst moments of the transfer eqnation, 

dtJ + d,W = r^^-XeJ + V^W[Xe + ed,Xe] 

~V‘^eder]e + ^V‘^e‘^dlr]e - \xeJV^ 

-^-K'^^V,V,e^dlxe 

and 

dtW + = -XsH^ + {2Ve - edeVe) + {Xe + ed,Xe) 

-^XsH^V^ - sd.Xe - \Q^^%Vue^dlxe • (B5) 

Blandford & Payne (1981a) added a photon sonrce term to the right side of their 
transfer eqnation, withont adding any corresponding absorption term. This is fnndamentally 
inconsistent with thermodynamics (if absorption is not inclnded, the photon source can 
never come into thermodynamic equilibrium with the radiation held). It is equivalent to 
including the terms on the right sides of equations ( |B3|) -(P5|) that involve the emission 
coefficient but neglecting all the terms that involve the absorption coefficient Xs- 

Neglecting the absorption terms compared to the emission term is a valid approximation 
only if the specihc intensity of the radiation held is negligible compared to the source 
function Se = Pe/Xe o-i M photon energies^ i.e., only if self absorption is never important at 
any energy, anywhere in the system. This is rarely the case in astrophysical systems. For 
example, it is not the case in the Comptonizing regions around accreting neutron star and 
black hole X-ray sources. 

Neglecting the absorption terms in the transfer equation leads to equations that have a 
diherent mathematical character from equations ( P3|) -( p5[) , because in this case the right 
sides of the transfer and moment equations do not depend explicitly on the radiation held. 


:JV^ + K^^ViVj 


edeXe 


(B4) 
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